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Abstract 

The Majorana representation of spin operators allows for efficient field-theoretical description of spin-spin correlation 
functions. Any N-point spin correlation function is equivalent to a 2N-point correlator of Majorana fermions. For 
a certain class of N-point spin correlation functions (including “auto” and “pair-wise” correlations) a further sim¬ 
plification is possible, as they can be reduced to N-point Majorana correlators. As a specific example we study the 
Bose-Kondo model. We develop a path-integral technique and obtain the spin relaxation rate from a saddle point 
solution of the theory. Furthermore, we show that fluctuations around the saddle point do not affect the correlation 
functions as long as the latter involve only a single spin projection. For illustration we calculate the 4-point spin 
correlation function corresponding to the noise of susceptibility. 
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Spin systems are notoriously difficult to describe using field-theoretic methods due to non-Abelian nature of the 
spin operators Often one tries to circumvent the problem by mapping the spins onto a system of either bosons 
or fermions, for which a standard field theory can be developed 0. There is no unique recipe for such an approach. 
Several formulations have been put forth for solving specific problems, including the Jordan-Wigner J3] and Holstein- 
Primakoff EQ transformations, the Martin a Majorana-fermion and Abrikosov 0 fermion representations, as well 
as the Schwinger-boson GHE) and slave-fermion mm techniques. 

The Jordan-Wigner transformation is the only “exact” mapping between the spin-1/2 and fermion operators as it 
preserves not only the operator algebra but also the dimensionality of the Hilbert space. However, it is a non-local 
transformation specific to one spatial dimension |2] where it is often applied to Bethe-Ansatz-solvable models or 
their variations (a generalization to higher dimensions does exist mm. but it lacks the simplicity of the original 
approach). All other mappings suffer from the following two problems: (i) the Hilbert space of the fermionic or 
bosonic operators (the “target” Hilbert space) is enlarged as compared to the original spin Hilbert space, and (ii) the 
resulting theory in the fermionic or bosonic representation needs to be treated perturbatively, which often leads to 
complicated vertex structures (see e.g. ED). The former issue may be resolved by additional constraints HUB ED or 
by projecting out unphysical states J6j at the expense of further complications such as the appearance of non-Abelian 
gauge fields ll9l flOlfl3lfl6l . 

The Majorana-fermion representation, suggested by Martin a, offers a possibility to avoid both types of prob¬ 
lems mentioned above: (i) The target Hilbert space is indeed enlarged, but merely contains two (or more) copies of 
the original physical spin Hilbert space mm. Matrix elements of physical quantities between different copies 
of the original Hilbert space vanish and thus the correlation functions may be evaluated directly in the target Hilbert 
space (this fact is often not fully appreciated; below we justify and further illustrate this statement), (ii) The Martin 
transformation 0 represents the spin operators in terms of bilinear combinations of Majorana fermions. The result¬ 
ing Majorana theory appears to be interacting (similarly to the situation with the Jordan-Wigner transformation 0 
and other fermionic representations EHOD) and /V-point spin correlation functions are equivalent to 2A-fermion 
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correlators. However, a certain class of /V-point spin correlation functions can be reduced to A'-point Majorana cor¬ 
relators ED El- In this case, complicated vertex structures do not arise. This significant simplification applies to 
correlation functions which involve only a single spin (“auto-correlation” functions) or an even number of spin op¬ 
erators for each spin (“pair-wise” correlations). In the present paper we generalize this approach to higher-order 
correlation functions and develop a technique for practical calculations based on the Keldysh formalism [23]. 

We illustrate our general arguments by a specific example. For a model of paramagnetic spins coupled to an Ohmic 
bath [24fl we determine the four-point correlation function describing the noise of susceptibility [25l . This quantity 
is closely related to the recently measured inductance fluctuations in SQUIDs ||26| . In addition, it is a direct measure 
of non-Gaussian fluctuations that are relevant also in other physical contexts EZM]. The noise of susceptibility is 
distinct from the better-known four-point correlator “second noise”, which always comprises a significant Gaussian 
contribution [3014321 . 

The paper is divided into two major parts. In Section[l]we introduce the Majorana representation of spin operators 
and address general issues related to this approach. In Subsection [L2] we discuss the fact that spin-spin correlation 
functions can be directly calculated as correlations of the Majorana fermions despite the enlargement of the fermion 
Hilbert space. In Subsection 1.3 we demonstrate the simplification of the theory for the case of auto- and pair- 
wise-correlations. In Subsections 1.4 and 1.5 we formulate of the above correlators within the Keldysh path-integral 
approach. 

In Section [2] we apply our general arguments to the problem of the noise of higher-order spin correlators in the 
Bose-Kondo model. Here we introduce the path-integral formalism in the Matsubara representation and show that 
certain higher-order spin correlators can be calculated in the saddle-point approximation. Further technical details are 


provided in the Appendices. In Appendix A we present the traditional diagrammatic perturbation theory. In Appendix 
[B] we justify the saddle-point approximation used in Section[2] Finally, in | Appendix C we analyze a remarkable gauge 
freedom in our model. 


1. Majorana Representation for Spin Operators 


1.1. The Martin transformation 

In this paper we focus on the following Majorana representation of the spin-1/2 operators introduced by Martin 
in 1959 0 in the framework of generalized classical dynamics: 


S ~ 2 e a:fiyt)l3Vy > 


s x = -ifjyfjz , s y = -ifj-f) x , S z = -ifjxfty 


The Majorana operators obey the Clifford algebra 


0 ) 


{f/a-rip) = , C =1 /2, (2) 

where {.,.} denotes the anticommutator. The Majorana representation has been used in a variety of physical contexts, 
cf. Refs. [2l I2T1 l22l l33l . (We remark that another normalization of the Majorana operators, if - 1, used in Ref. [211 
only changes some numerical prefactors at intermediate stages of the calculation.) 

The above representation with the real Majorana fermion operators, 

ff = (3) 

perfectly reproduces the SU(2) algebra of the operators S a 

[S a ,S' 3 ] = ie a p y 8*, (4) 

and explicitly preserves the spin-rotation symmetry. 

Applying standard field-theoretical methods to fermionic systems implies the existence of a Fock space. In a 
faithful representation of a spin 1 /2, the dimensionality of the fermionic Fock space should coincide with the di¬ 
mensionality of the Hilbert space of the spin. For example, the Jordan-Wigner transformation represents a system 
of N spins-1/2 with the 2 v -dimeiisional Hilbert space in terms of N fermions with the 2 JV -dimensional Fock space. 
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However, each Jordan-Wigner fermion may be expressed in terms of two Majorana fermions. In contrast, the Martin 
transformation (|T|) represents each spin in terms of three Majorana fermions and hence does not preserve the dimen¬ 
sionality of the spin Hilbert space. 

One method of dealing with this issue is to express the Majorana (or “real fermion’') 77 -operators in terms of 
’’complex” or ’’Dirac” fermions (this is a common slang used to distinguish usual fermions from Majoranas; these 
fermions do not necessarily obey the Dirac equation) and use their respective standard Fock spaces. This requires an 
even number of Majorana fermions. The simplest possibility is to add one auxiliary Majorana (cf. the drone-fermion 
representations of Refs. 151 [20113411351 ). In this case, the resulting target Hilbert space is 4-dimensional and can be 
interpreted as two copies of the spin Hilbert space GOlETl . To make the spin-space isotropy even more explicit, one 
could add (and pair) an extra Majorana to each of the three t) a \ this, however, results in an 8 -dimensional Hilbert 
space, equivalent to four copies of the spin Hilbert space. In the following subsection we demonstrate that regardless 
of the construction of the Hilbert space, spin correlation functions coincide with Majorana-fermion correlators. 

1.2. Equivalence of spin and Majorana correlation functions 

Now we show explicitly that correlation functions of spin-1 /2 operators can be directly computed in the Majorana 
representation 0 (regardless of the explicit construction of the Majorana Hilbert space). Indeed, the equation ([TJ 
gives a representation of SU(2) with S = 3/4. Thus, it is necessarily a direct sum of a certain number of irreducible, 
two-dimensional spin- 1/2 representations, i.e, we obtain an integer number of copies of the spin (the actual number 
is determined by the particular choice of the number of auxiliary Majorana operators). The spin operators (jTJ do not 
switch between the copies, thus any trace of an operator built out of spin operators 0 is given by the number of copies 
times the trace in a two-dimensional spin space. This leads directly to the desired result as we discuss in more detail 
below. 

Let us remind the reader that any function of the spin-1/2 operators (i.e., of Pauli matrices) is in fact a linear 
function. This follows from the following relation for the spin-1/2 operators 

S a S fi = l -e afiy sy + ] -6 a p. (5) 

For an arbitrary number of operators describing the same spin this implies that 

S a '.. .S a " = a 0 S° + ao, ao,a x ,a y ,a z e C, (6) 

where the complex numbers a (h a x , a y , and a z depend on the sequence {a,}. Since spin operators are traceless, the 
trace over the 2-dimensional ids = 2) spin Hilbert space Tn, of the above combination of spin operators is given by 
the constant term oq: 

Tr s {s ffl ...S a "} =d s ao(a u ...,a„), d s = 2. (7) 

The Martin representation (|T]> preserves the commutation relations as well as the relation <|5j, and hence the 
equality (| 6 | remains valid in the Majorana representation (where the spin operators should be understood as Majorana 
bilinears). In particular, the coefficients a, remain the same. The trace over the Majorana Hilbert space can be 
performed by noting that Majorana bilinears are traceless due to their anticommutation properties. Again [cf. Eq. (jTJ], 
the only remaining contribution is given by the constant term a q: 

Tr M {5“>..J a "} = d M a 0 (a u ...,a n ), ( 8 ) 

where is the dimension of the Majorana Hilbert space. Thus, tracing an arbitrary product of the spin-l/2 oper¬ 
ators over the spin and Majorana Hilbert spaces yields the same result up to a numerical factor, determined by the 
dimensionalities of the Hilbert spaces). 

Similar arguments were put forth in Ref. Il20l in the context of the drone-fermion representation. The results of 
this Subsection were implied in Refs. f2Tl[22ll and given without proof in Ref. |[25l . 

The above statement can be readily generalized to an arbitrary ensemble of spins. Indeed, any function of spin 
operators is still linear in the components of each spin. If the operators on the left-hand side of Eq. 0 describe more 
than one spin, then on the right-hand side additional terms appear, which contain all possible products of operators 


3 


related to different spins (for example, in the case of two spins the right-hand side reads « ( > + a„S + b„S “ + c a pSfS^). 
However, all such additional terms are still traceless, and hence the only change in Eqs. 0 and ([8]) will be in the 
constants ds and dM- 

Consider now a real-time spin auto-correlation function: 


(§ a '( tl )... 5 “"(?„)) 


Ti m 


(9) 


where p = exp (~/3H) is the non-normalized Gibbs density matrix. If in addition the spins are coupled to other degrees 
of freedom denoted collectively by X (this is also described by H), then following the line of arguments presented 
above we can write p = po(X) + S a p a (X), where po(X) and p a (X) are matrices in the X space. Moreover, po(X) is 
the reduced density matrix describing the rest of the system (the X-degrees of freedom). In this case, the partition 
function may be written as 


Tr |p) = Tr x Tr s (p) = Tr x Tr s [p Q (X) + S a p a (X)) = d s Tx x (p 0 (X)} . (10) 


Now, each Heisenberg spin operator S a ‘ (ti) is related to the Schrodinger operators (i.e., the Pauli matrices) by the time 
evolution operator U(t,t') = exp[-f J r drH(T)]. The latter can be expanded similarly to the density matrix. Thus we 
can write 5 Ql (fi)... S a,, (t n )p = Ao(2Q + S^A^X), where Aq(X) and A a (X) are matrices in X-space. As a result, we 
can formally perform the trace over the spin variables in the auto-correlation function (|9]i and find 


(s a '(t l )...S a "(t„)) s 


TrxTrs {a 0 (X) + 5“A ff (X)} 
^sTr x IpoW) 


Tr x (a 0 (X)} 
Tr x (poW) 


( 11 ) 


In the Majorana representation the structure of the above equations remains the same. The averages can be 
formally calculated similarly to Eq. 0- As a result, we arrive at the expression 




Tr x Tr M {a 0 (X) + S n A a (X)} 
^AfTr x {p 0 (X)} 


Tr x {a 0 (X)} 
TrxfpoW)’ 


( 12 ) 


which is identical to Eq. CD- The generalization to the case of general (multi-spin) correlation function is straight¬ 
forward (see above). Thus we have demonstrated that spin correlation functions can be calculated with the help of the 
Majorana representation 0 without any projection onto ’’physical” states: 

(s ai (h ).. J ff "(f„)) s - (s ai (t 0. -.S a "(t n )) M . (13) 

The operators on the right-hand of Eq. represent the Majorana bilinears (JTJ. 


1.3. Simplified representation for autocorrelation functions 

The arguments presented in the previous Section show that any /V-point spin correlation function may be repre¬ 
sented in terms of a 2A-point correlation function of Majorana fermions. Here we demonstrate that for a certain class 
of correlation functions this correspondence can be significantly simplified [ 2J_. 22 36J. 

Let us rewrite Eq. 0 as follows 

S° - @f/ a , 0 = -2i?i x f) y fj z , Q 2 = 1/2. (14) 

One can easily verify that the operator 0 commutes with all three Majorana operators ij a . Consequently, this operator 
also commutes with any Hamiltonian expressed in terms of ?/,,, and thus the corresponding Heisenberg operator is 
time-independent. The averaged product of a pair of spin operators can now be represented as follows 

(s°(o)SV)) M = (®ut)®wo) = \ {umt')) ■ as) 
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Thereby a two-point spin correlation function reduces to a two-point (rather than four-point) Majorana-fermion cor¬ 
relation. 

Unfortunately, the above relation cannot be directly generalized to time-ordered correlators (or Green’s functions) 
due to the fact that spin and Majorana operators are influenced by time ordering in different ways. Explicitly, a 
time-ordered average of two spin operators is given by 



(s a (t)§e(o), 
(§P(f )§*(!)), 


t > t' 
t < t' 


1 , t>t' 

2 \{fl/i(t')f] r y(t)) , t < t’ ’ 


( 16 ) 


where the latter expression differs from the time-ordered average of the two Majorana-fermion operators by the ab¬ 
sence of the minus sign in the lower line. 

While this problem can be circumvented by introducing Green’s functions of the operator 0 If37l . we consider here 
a simpler approach. The missing sign can be compensated for with the help of an auxiliary Majorana fermion m that 
anti-commutes with the three operators fj a , i.e., = 0 and m 2 = 1/2. Since any Hamiltonian will be expressed 

in terms of bilinears of fj a , the operator m commutes with the Hamiltonian and, thus, is time independent. We keep, 
however, its formal time argument in order to be able to treat time-ordered operator products correctly. We notice that 


( T 


{im(t)r] a (t)im(t')r]i j (/ 



{im(t)rj a (t)im(f )Tj P (t ', 


t > f _ 1 f( TJ a (t)Tlfi(t ')) , t > t' 
t<f 2 \{rip{f)T]a{t )), t<f 


Thus we arrive at the identity 


(t {s°(osV))) m = {t \immxt) mnuo}) ■ 


(17) 


Here the 2-point spin correlator is again expressed in terms of a 4-point Majorana correlation function. However, in 
contrast to the direct application of the Martin transformation 0- here two of the Majorana operators m(t) and m(t') 
do not have any dynamical properties and only serve the purpose of writing the time ordering (16 1 in a compact form. 


Similarly, we can use the auxiliary operator m(t) to express higher-order correlation functions. For a 4-point spin 
correlator we find 


(t {s“(f I )s^(f;)5 5 fe)s’'(4)}) M = (t {muoUhm', )UA )«(f 2 )7fe(? 2 )«(4)7r(4)}) • 


( 18 ) 


The correlation functions ( fT7| and ( [T8| are in fact auto-correlation functions in the sense that they involve opera¬ 
tors describing the same spin. Clearly, the simplification (15 1 cannot be extended to different spins since the operators 
0i and 02 anti-commute and 0!02 4- 1. Therefore, at the level of 2-point correlation functions the simplification 
described in this section applies to auto-correlation functions only. Generalizing this technique to higher-order corre¬ 
lation functions, we can compute autocorrelators, such as the 4-point function ( fl8] >, as well as “pair-wise” correlators 
comprised of pairs of operators for each spin, such as 

(7’{5^(f 1 )5f(f;)5^ 2 )5^')}) M . 

Other correlators, such as (t {5“(f)5^(f')})> have to be computed by different methods. 


1.4. Spin correlation functions in the Keldysh formalism 

Real-time correlation functions at finite temperatures can be conveniently computed within the Keldysh formalism 
m EH. The calculation amounts to finding the generating functional Z A j23l and then taking the derivative with respect 
to the source fields. 

For a spin system, the generating functional may be defined as follows 

£>[...] exp jiS 0 + ijdt (a^SI + , (19) 
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where D[... ] denotes the appropriate measure of integration whereas So represents the action of the model under 
consideration. In particular one could choose to integrate over the SU(2) group manifold and D[... ] would represent 
then the appropriate Haar measure (T). I' 1 this paper we choose a more straightforward method of integrating over 
real Grassmann variables representing the Majorana operators of (JT]». In ( 19 1 A c ^ q) are the source fields. The super¬ 
scripts cl and q refer to the “classical” and “quantum” variables (23) that are defined as the sum and difference of the 
corresponding fields belonging to the upper ( u ) and lower (d) branch of the Keldysh contour 

s c ™ = ^=(s“±s d a ), 4 ' w = ^( t “ ± 4 ). 

The ’’classical” source term defined in this way describes the physical probing field, A C J = y2B a , while the ’’quantum” 
term is only needed to construct the correlation function and is set to zero at the end of the calculation. 


Taking the derivative of the functional (jl9b with respect to the source fields A c J (q \ one finds the spin correlation 


functions. In particular, the one-point function defines the magnetization 

<)Z 


V2AT = ( S a ’ c, (t )) = -i 

The spin susceptibility is given by a 2-point function 


SA^t) 


XadSt') = i(T K S a ’ d (t)S^(t f )) = -i 


.1-7=0 


b 2 Z.i 


8AP’ cl {f)6A a *(t) 


( 20 ) 


( 21 ) 


i«=0 

while the noise spectrum ED is determined by a different 2-point correlator, (TKS c J(t)S c J(t')y 

Below in Section [2] we will be interested in a specific 4-point function that determines the experimentally acces¬ 
sible noise of susceptibility |[25l(26) 


C x {t u t\,t 2 J' 2 ) = 

and we focus on the case of identical spin indices. 


b 4 Z.i 


bd“- c, (f;)b/i Q '-9(f 2 )b/i“- w (f;)bT Q ''«(f 1 ) 


( 22 ) 


i-7=0 


1.5. Correlation functions of Majorana fermions in the path-integral representation 

Let us now reformulate the mapping between the spin and Majorana fermion operators discussed above in Sections 


1.2 and 1.3 in the language of the Keldysh functional integrals. In the operator language we have established the 
correspondence between the spin /V-point functions and Majorana 2Appoint functions. In special cases (namely, 
for autocorrelation functions and ‘pair’ correlators) we found simpler relations 0 and ( p~8j ). 

For usual (Dirac) fermions the path integral approach is based on the concept of coherent states [DIIII, which 
are eigenstates of the fermionic annihilation operators. An annihilation operator can only be constructed from two 
Majorana fermions. In many-body problems, this is typically achieved by either “halving” Il38lj40l or “doubling” 
eseded the number of Majorana operators. Keeping in mind applications to systems with a small number of 
degrees of freedom (such as the single-spin problem discussed below), we adopt the doubling procedure, where we 
add a free Majorana fermion to each of the operators introduced by the Martin transformation 0- The interaction part 
of the Keldysh action can then be formulated in terms of the Grassmann variables q„ corresponding to the Majorana 
operator q a in the previous sections. The additional Grassmann variables appear only in the “free” quadratic part of 
the action and are decoupled from the physical system under consideration. 

The Majorana-fermion 2/;-point correlation function on the right-hand side of Eq. ( p~3j ), or rather its time-ordered 
counterpart, can be obtained either by 2 n differentiations of the generating functional with respect to Grassmann source 
fields each coupled to a single Grassmann variable q a , or by n differentiations with respect to c-number source fields 
coupled to pairs of Grassmann variables q„. These pairs should then be chosen to represent the spin components 
according to Eq. 0 as in For the autocorrelation functions or the pair correlators we would like to use the 


simplified correspondence, i.e., Eqs. (17 i and (18 i. Here the auxiliary Majorana fermion m can be represented either 


by one of the free Grassmann variables used to construct the functional integral, or by one out of yet another pair of 


Grassmann variables that are added to the theory specifically for the purpose of computing the correlators (17 1 and 
( fi~8j ). The source fields can again be chosen as either Grassmann variables or c-numbers. In the explicit calculation 
below we chose the latter option. Of course, physical results are independent of these technical details. 
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2. Spin correlators in the Bose-Kondo model 

In this section we apply the general conclusion reached in the previous section to a specific example. For simplicity 
of the presentation we use, first, the Matsubara technique, while the final results are formulated in real time in the 
frame of the Keldysh formalism. 

As a model we choose a zero-field spin-isotropic Bose-Kondo model (see Ref. Il42l |43) and references therein). 
This model appears in various physical contexts, including spin glasses and liquids [S 44l [45l . In the context of 1// 
noise a similar model is known as the Dutta-Horn model Il24ll . where a large number of independent spins are coupled 
each to their own bath. In the Majorana representation introduced above Q. the model Hamiltonian reads 

H = $X + H b = - l -X a e a , iyWl7 + H b , (23) 

where H B is the Hamiltonian of the bosonic bath controlling the free dynamics of The latter may be characterized 
by a Matsubara correlation function {TX a {T)Xp(r')) = 6 a pYl(r - r'). The function n(r - r') can be written as 

A 

r dx p( |x|)signx 

n = -T- , (24) 

J n x — iu> m 

-A 

where p(|x|) is the bath spectral density and a> m = 2nmT. In the Ohmic case considered here, p(\x\) = g\x\. For 
u> m <s A, this gives I l(ioj m ) ~ (2g/n)A - g\u> m \. One can perform an RG procedure by integrating out energies of 
the bath between A/b and A. As a result the coupling constant g is rescaled. The RG differential equation reads 
dg/dlnb = -2g 2 /n (see, e.g.. Ref. El ED). One can supplement this with the scaling equation for the quasi¬ 
particle weight d\nZ/d\nb = -2g/n, which could be important |461 as we are interested in Green’s functions of the 
Majorana fermionic operators. For go = g(ln b - 0) <s: 1, the renormalization effects are not important as long the the 
temperature is high T » Tk = Aexp[-7r/(2go)]- Here we assume this to be the case. Thus we can safely reduce the 
cutoff A (b) to a value of the order of temperature. 

2.1. Matsubara path integral 

We use the Matsubara imaginary-time technique (t = — ir, d T = -id,). The partition function (Z = f D[... ] exp [iS ], 
for brevity we omit the source fields) reads 


Z = 


l/T 

J D[$]D[r, a ] exp < / S B + i J" dr 


^ Ua(T) i d T ij a {j) + ^ X a (T)e a p Y Ti fi (T)ri 7 (T) 


the matrix M is of the form 


M = 


0 1 G 
1 0 1 
1 1 0 


Next we decouple the quartic Majorana-interaction in a different channel. To this end we rearrange 
is M.mtldn] = ~ J' drdr' M a pYl(T - t) rj a {T)T]p(r) ^(t'^/jCt') 

J " dTdr' n(r - r') [rj a (T)r] a (T')\ M afj f^W^/r')] . 


1 

4 


(25) 


Here S B is the free bosonic action. The first step is to average over the fluctuations of It, yielding 

Z = JD[r] a ] exp j-i j r dTTja(T) d T rj a (T) - 1 J'drdr'M^ n(r - r') rj a (T)rjp(T) /7„(t')?7/?(t')| , 


(26) 


(27) 


(28) 
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We now employ the Hubbard-Stratonovich transformation by introducing fields X„. These fields inherit the symmetry 
of Majorana propagators, therefore Xq,(t, r') = T ). The new effective action reads 




■ j' cItcIt' 


1 , , 1 £ ff (r, t') M 1 X^r, t') 


The Majorana Green’s function in ([29]) is 

(G-') tt , =S(T-T')d T ,-Y a (T,T'). 

The function 1 T(t - t') is positive and non-zero. The standard form reads 

gxT 2 


rf(r - /) = 


sin"(7rr(r - r')) 


(29) 


(30) 


(31) 


It is cut off at short times, |t - r'| < l/A, leading to the maximal value of order gA 2 . We use the divergent form 
keeping the regularization and renormalization in mind. Diagonalizing M -1 , we choose the eigenmodes with positive 
eigenvalues to be real and eigenmodes with negative eigenvalues to be imaginary, such that the overall sign of the 
action term in the exponent is negative, and the functional integral over X„ converges. In other words, we choose 
X = X' + zX", where the real part is ‘diagonal’, X' = X' • (1,1,1), and the imaginary part X" is orthogonal to it, 
(i.e., X" + X" + X" =0 and X" describes two degenerate modes); with this choice the 3D integral over X' and X" 
converges. The redecoupled action ( [29] ) is again quadratic in Majorana Grassmann variables q, r , which allows us to 
integrate them out and to obtain an effective action of X-fields: 


iS [X„] 


k Z Trl08 ( G "')'i/ 


tirdf 


a={x,y,z] 


, £ ff (r, t 0 Mj X^t, t') 

n(r-r') 


(32) 


2.2. Saddle point solution 

We can now identify the saddle point and fluctuations of the effective X-action. The saddle-point solution is found 
by expanding X ff = X( )r , + S£, r . The linear order in <)X vanishes for 


£oq-(t - t') = n(r - T r )M a pG(^(T - t) , 


(33) 


where Go/? = G/?[Xo/?L A straightforward calculation, in which we disregard the broadening of Go/? on the r.h.s. of 
Eq. ( [33] ), leads to 

A 


^0 orO^n) — 


f 


dx gxcoth[jSx/2] 
n x- ie n 


-A 

Upon the analytic continuation ie n —> e + iO and for e —> 0, we obtain the retarded self-energy 


(34) 


X£,(e -> 0) = -2 igT = —iF . 


(35) 


Here we recognize T = 2gT to be the (Korringa) relaxation rate. In terms of the eigenmodes this solution means 
Xq = 0, whereas X£,(ie„) = X 0 „(iT„). 


2.3. Fluctuations 

To study the fluctuations we expand the trace-log-term in ([32]) to second order in dX. The action reads 


iS as 


1 

4 



dr i dr 2 dr 3 dr 4 G 0 , a (T\ - r 2 ) dX a (r 2 , T 3 )Go, a (T 3 - r 4 ) 5 Xq.(t 4 , n) 


1 r tiZ a (T, T ) M * SI.fi(T, T ) 

- I dr dr -—--- 

4 J n(r - t') 


(36) 
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The Fourier transform of 61 is introduced via 

61 a (e„,v m ) = J dTdT'e iv '^ T+T,) e ie ^- T,) 6Y a (T, /). 


( 37 ) 


The direct analysis shows that one of e„ and v m must be fermionic and the other bosonic, so that both e„ + v m and 

n 

T 2 


iS 


SL 


T 3 

T 


ZZ G(),a(.€n Vm) ^0,a^n V m ) S^ a (^6 n ,V m ) S^ a (^€ n , Vm) 

a e n ,v m 

2 XI A ( £l + €2) V ) M aP -V) . 


a (3 6162 V 


Here 


A(ia> m ) = dr 


1 IT 


1 


n(r) gnT 2 


l/T 

f 


dr e" 


sin 2 ( 71 T r) = 


2g7tT 2 


hn, 0 


2m,1 


2m,—1 


(38) 


(39) 


The mean-field Green function reads 


Go, a (ie„) — 


(40) 


Since e„ is necessarily fermionic, we have |Go, ff 
second one which is proportional to 1 /g. This important observation allows us to disregard the first term of 


1 

ie„ + /T sign(e„) ' 

1/(7tT). Thus, the first term of (38 1 cannot compete with the 

and 


essentially all the contributions of the second and higher orders, originating from the trace-log term of (32 1 . This in 
turn simplifies calculations of the higher-order spin correlators in the next section. 


The argument above for the smallness of the first term of (38 1 is based on the discreteness of the Matsubara 


fermionic frequencies. Ultimately, we are interested in real times and the behavior in various frequency ranges, 
including low frequencies u> <K T; hence one should be careful with the estimates. The expressions above indicate 
that for e —» 0 the Green functions Gy ,, in the first term of ( |38|> might become of order 1/F = 1 /(2gT). This, in turn, 
might imply that the (prefactor of 61? in the) first term of (|38[) scales with g 2 and dominates over the second term 


g . To clarify the situation in the low-frequency range we perform a direct Keldysh calculation in Appendix B 


We conclude that the first term of (38 1 does not become large in the domain of low real frequencies. Thus, the second 
term of (f38l) dominates. 


2.4. Averaging over fluctuations 

The knowledge of the propagator of the r)X-f1uctuations allows one to construct a new perturbative series, starting 
at the fixed-point solution given in Eq. ( |35j ). The new perturbative expansion for the Green function is based on the 
following series: 

G a (T,fl) = G 0 ,ff( t, t') + J' dTidT 2 Go,a(T,Ti)61 a (TuT 2 )Go,a(T2,T') 

+ j" dTidT2dT 3 dT 4 GoA T ’ Tl)6t a (TuT 2 )Go,a(T2, T 3 )dt a (T 3 , T 4 )Go, q .(t 4 , t') + ... (41) 

This series is to be averaged over the fluctuations to obtain (G a (t, 0) 5Y - However, on the right-hand side the spin 
index a is the same in all the terms, essentially because the saddle-point Green function is diagonal in spin space. 
Taking into account the dominance of the second term in the action ( [38) , we conclude that 

(42) 


(61 a (...)61 a (...))^M aa = 0. 


because the diagonal entries of the matrix M are zeros (note that in Eq. ( |42) the contributions of the real and imaginary 
E-components cancel each other, cf. the discussion above Eq. ([32)). Thus, the fluctuation-averaged Green function 
coincides with the saddle-point solution 


( Ga(t , = Go,ff( f > t '). 


(43) 
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Figure 1: The diagram for spin susceptibility. 
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Figure 2: The ’lowest order’ diagrams to non-equilibrium noise of susceptibility. 


Moreover, for averages of Green-function products with the same spin index a we can substitute all G a with their 
saddle-point values due to the relation ( |42| >; this simplifies calculations considerably. In the next section we employ 
this convenient property. 

2.5. Correlation functions 

As a first example, we calculate the diagonal spin susceptibility 

X(t,f) = i{T K S a ’ c, (t)S a ’Hf )) = -i(t k {m(t)fla(t)] cl . (44) 

The only contribution to x(t, f) is given by the diagram in Fig. [T] We immediately arrive at the standard result 


X M = ~ 


1 T tanh w<<7 - 1 cjT 


2 w 2 +r 2 


4 T to 2 + r 2 ' 


(45) 


Notice that no vertex corrections appear in our case, and the precision of this calculation rests solely on the precision, 
with which the self-energy is evaluated. 

We further demonstrate the power of our approach by calculating one of the higher-order spin correlators. We 
evaluate the connected part of the 4-th order correlator (|22|, which is related to noise of spin susceptibility 


Using (18 i we obtain 

C x (t A j\,hJ 2 ) = -(7>S“- c, (f 1 )S^(f , 1 )S“- c, (f 2 )S a - , « , 2)) 

= - (T r: [m(t\)fi a (,h)] cl \m{t\)T] a (t\ )\ q [m{t 2 )f) a (t 2 )] cl \ih(f 2 )r] a (t' 2 )]‘ l 'j . 


(46) 


The discussion in the previous section implies that the connected part of this correlator is given by the six diagrams 
depicted in Fig. |5] Here, the double solid lines stand for the saddle-point Green functions Go,,,, whereas the dashed 
lines represent the trivial correlators of the conserved quantity m. As shown in the previous section, we may use the 
saddle-point Green functions, since the contribution of the d£„-fiuctuations vanishes, cf. (42 1 . For completeness, we 
show the result, which is put into context and interpreted in detail elsewhere 


C c y om (v, 0 J U 0 J 2 ) = 


iT 2 


u> i + to 2 + 2zT 


8 T 2 (oj [ + i\ )((o 2 + /I )(oo\ + v + iY){oj 2 — v + /I j 


If a spin correlation function involves different spin components, one can no longer rely on the saddle-point 
contributions. One example is the correlator 


(r K s a ’ c \t i )s a ’ q (t\ )s^ c \t 2 )s^\t ' 2 )), 


(47) 


related to the correlations of susceptibilities in different directions. The non-zero off-diagonal fluctuations ^ dt.adl.pj 
contribute to this correlator, and thus additional diagrams have to be considered. This is, however, beyond the scope 
of the present paper. 
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Figure A.3: The leading contribution to the self-energy. 


3. Conclusion 

In this paper we have demonstrated the efficiency of the Majorana representation of the spin-1 /2 operators <|T]». We 
have shown that the representation 0 allows for a direct calculation of spin correlation functions in terms of Majorana 
fermions despite the fact that their Hilbert space is enlarged as compared to that of the spins. The precise construction 
of the Majorana Hilbert space was shown to be irrelevant as far as the correlation functions are concerned. 

The Majorana representation is particularly efficient in the case of auto-correlation functions (or “pair-wise” cor¬ 
relation functions). Such /V-spin functions can be represented in terms of /V-point Majorana correlators, which signif¬ 
icantly simplifies calculations. In particular complicated vertex structures do not appear. 

As an example we have revisited the well known Bose-Kondo model. We have developed the Keldysh path-integral 
approach and have shown that the spin relaxation and susceptibility are efficiently described within the saddle-point 
approximation. Moreover we have shown that correlation functions containing a single spin projection can also be 
efficiently calculated at the saddle-point. In particular we have evaluated a 4-spin correlation function corresponding 
to the noise of susceptibility. 

It would be interesting to apply our approach to a wider range of physical problems, for example, to a sub-Ohmic 
Bose-Kondo model describing physics of spin glasses Il42ij45l . 
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Appendix A. Perturbation theory 

Given the Majorana representation for spins <|TJ, the identities ( fT7| , ( fj~8j >, and the Hamiltonian ( f23] >, one could 
calculate correlation functions within perturbation theory by diagrammatic expansion. Here the zeroth-order Hamil¬ 
tonian consists of the bath part Hg only. The dressed Majorana Green function is self-consistently constructed, using 
the lowest-order self-energy. The lowest-order self-energy diagram is depicted in Fig. A.3 it contains the “free” 
Majorana Green function G/, ff (<A), 


G**(to) = (fi}±iOT\ (A.l) 

K(OJ) = - f ^(n K (oj + Q)GjjQ) + n Jf (aj + Q)Gf a (Q)) = -2igT^l+0^, (A.2) 

G R (u) = (oj- 2*(w)) _1 = (to + dT 1 , F = 2gT. (A.3) 

Here the Keldysh components of the bath correlation functions read 

n %(f, f) = <v (' r K x a ' a (t)x a - h (f )) = n ab (t, f), (A.4) 

n r '\cj) = ± g oj - iD, n*( w ) = coth(^) (nV) - n A ( W )) , (A.5) 


li 





Figure A.4: Examples of diagrams, which contain connected Majorana propagators. 



Figure A.5: Dyson equation for the partially dressed bosonic interaction. Thick solid lines are dressed Majorana propagators, wavy lines correspond 
to the bosonic bath. 


where D = (2g/n)A, the spin indices a,/3 = {x,y,z}, the Keldysh indices a,b = {cl, q], and the ’’classical” and 
’’quantum” operators X a ’ cl(?) = (X aM + X a - d )/ V2. 

Typically, this method is used to calculate 1- or 2-point correlation functions, e.g., magnetization or susceptibility. 
In lower orders Majorana propagators are never connected by a bosonic line. Such an element only appears in higher 
orders or in higher correlation functions as discussed in (25). An example is shown in Fig. A.4 Here it is in principle 
necessary to consider more complex diagrams including, e.g., ladders of bosonic propagators. In the perturbation 
theory with a vanishing unperturbed spin Hamiltonian it is by no means justified to simply neglect these diagrams. 
As an illustration let us consider renormalization of the interaction line. The dressed interaction carries four times (or 
frequencies), the Keldysh indices a,b,c,d = {cl,q}, and the spin indices a,[3,y, 6 = {x,y,z}', diagrammatically it is 
depicted by 


Y<ab,cd 
^ a(5,y8 


a, a 


c, 7 



(A.6) 


For demonstration, we pick one of several possible components of f that carries identical spin indices on the left (as 
well as on the right) and the Keldysh indices (21,21), which correspond to the retarded-retard ed co mponent, that is 
T RR . The partially dressed interaction is obtained from the Dyson-type equation, depicted in Fig. A.5 by summation of 
the leading contributions in the small-g expansion. These are constructed out of the bosonic correlator 1 \ K ~ 2 gT = I 
combined with Green’s functions G R G A such that the upper and lower halves of the complex plane each contain one 
of the Green-functions poles. These leading contributions are of the same order in g as the bare bosonic line, thus we 
might suspect a strong renormalization of the interaction line. This is, however, not the case due to a cancellation. 
Taking into account both contributions depicted in Fig. A.5 we obtain 


r y) = n K (v)M aP n V - Q)M ay G R (co 1 + G)G *(« 2 - G)r rr (<v u cj 2 , O) 

+ f ^ n7v - O )M ay G R (to l + O )G r (cj 2 - Q.)r™(co U 0 J 2 , -n - + 0 J 2 ) (A.7) 

The Green functions are calculated within the approximation I. R (a>) = -2igT = -;T, which is justified and consistent 
in the high-temperature regime T » T K , v, oj\ , w 2 , T. Let us assume that the leading term in 1 A>A ' does not depend on 
the third frequency and splits T RR = I + F RR . If this is the case, the integrals in the second and third term become 
equal, and we find that the renormalized interaction coincides with the bare one: 

= Tl K (0)M a p. (A.8) 
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The problem in this perturbative approach is that the non-perturbed spin Hamiltonian is zero (the non-perturbed 
Hamiltonian consists of the bath Hamiltonian Hb only). In other words, the spin-bath interaction is not weak as 
compared to the energy scale of the unperturbed spin dynamics. One might therefore expect that higher-order diagrams 
are important. Having observed a particular cancellation we cannot in principle exclude other important higher-order 
contributions. Instead of evaluating multiple higher-order diagrams we choose the path-integral approach in the 
following section allowing for a more straightforward analysis. 


Appendix B. Majorana path integral in the Keldysh representation 


In Section [2] we have developed a path-integral technique, which allowed us to obtain the self-energy (Korringa 
relaxation rate) as a saddle-point solution. We have argued that the fluctuations around this saddle point can sometimes 
be disregarded. Namely, this is the case for a correlation function involving only one spin component. This conclusion 


was based on the smallness of the higher-than-linear contributions to the trace-log term of the action < 32 1 . We have 


shown this using the small parameter T/e„ for the fermionic Matsubara frequencies. To be able to treat also low real 
frequencies, we provide here the Keldysh version of the path-integral calculation. 

For the model considered here, m the partition function reads 

Z = J' D[X]D[rj a \ exp i^iS B + J" dt (rj a (f)id t T] a {f) + iX a {t)e a ^p p {t)p y (f))^ 

= f D[f ]D[rj a ] expjiS B + £ f dt (r, a (t)i^d t r, a (t) + } , (B.l) 

^ l a-Lul O-co I 


where the Keldysh index a takes the value u at the forward part of the contour and d on the backward part of the 
contour. At the first step we average over the fluctuations of X, which yields 

Z = f exp f . (B.2) 

Here (TX‘‘(t)X b (f)) = d n f : U ah (t - t'), and a , b are the Keldysh indices over which summation is implied. As in 
Section[2] we decouple the quartic term in a different channel 

iS M ,int[ria] = ^ J dtdt' M a fjTj U ah (t, t')T b rf tja if tfi rj b , a r] b 

= ~lf dtdt ' T i f ') A ^llaAjy) (B.3) 


Appendix B. 1. Qualitative considerations 

Prior to performing the full fledged Keldysh analysis, we provide here a qualitative argument based on the locality 
of the bath correlation function n ah (t - t') on the relevant time scale of order 1 /T. On this time-scale we can safely 
replace all components of Yl ab (t - t') by its classical (Keldysh) part, i.e., n flfe (f - f) ~ TT(?| - t 2 ) = (l/2)FL K (t - t'). 


This allows us to proceed similarly to the treatment in the Matsubara case in the main text (cf. Eq. 32 1 , and we obtain 


iS[Z a ] = - 


Z Tri °s( G « 1 )-^/ 


- I dtdt' 


?? b (tj')T a i T b M-l'Lf{t,f) 


a={x,y,z} 


n (t-n 


(B.4) 


where 


(c; 


-1 \ab 

, )tr = lT i W ' 


t')d,,-x ab (.t,t'). 


One can find the saddle point and again obtain the relaxation rate T = 2gT. This is done below in the full Keldysh 
calculation. Here we concentrate on the fluctuations 51L ab . On the relevant time scales (~ T _1 ) the function ft is local, 
ft(r — t') ~ 2gT6(t - t'). (Note that the delta-function should be understood as such only at relatively long time scales. 
For instance, it does not force us to take the Grassmann variables in (B.3 1 at coinciding times.) This locality means, in 
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turn, that the fluctuating self-energies 6T. ah are local. On the other hand they are anti-symmetric, "L ah (l, t') — —T, ba (t', t). 
This strong constraint implies that only off-diagonal (in Keldysh indices) components of £E fluctuate: 6Y. ud (t, t’) = 
-6Y. du (t',t) ~ 6(t - t'), whereas bE"" = 6I. dd = 0. Upon the Keldysh rotation (see below), this means that only 
the retarded and advanced components of 61. fluctuate. In addition the Keldysh component of the Majorana Green 
function G* can be neglected, since it scales as oc tanh(w/27')d(w). As a result, upon expansion of the trace-log term 


of (|B.4|), the only terms that can appear are of the type Tr[G R bE R G R c5E R G R (5E K ... ] and Tr[G A bE A G A £E A G^E A ... ]. 


Since 6Y. K ' A are local in time, these terms vanish. Thus we are allowed to disregard the trace-log term of (B.4 1 . Notice, 
that this argument is not valid at finite magnetic fields (B > T, see below). 

Appendix B.2. Full Keldysh calculation 


We now go back to the full Keldysh version of (B.4 1 keeping all Keldysh components of IT' & . One can decouple the 
quartic Majorana interaction with the help of complex bosonic fields Q„ via a Hubbard-Stratonovich transformation. 
The fields Q a inherit the symmetry of the Majorana propagators, therefore Q“ h (t, t’) = -Q br ,‘(t', t). Since the correlator 
Yl ab (t - t ') may have a complicated time-dependent structure, we choose to keep it in the numerator: 


is [t]a, Q a ] = Jdtdf y-plit) (G~X tfjf ) - ^Tf^n ab (t, f) Q ab {t, t') M~lQ ab {t , f) 

The Majorana Green function in ( | B .5 [ ) reads 

{G^)° b = iT ab 6{i-t')d t ,-v; b {ut'), 

After the standard Keldysh rotation, 

G = LGL = 


(g k 

g r \ 

A In* 

n R \ 

(g a 

o • 

n-Lnu- 

o) 


Z = LZL = 


0 


E a 

y. k 


“f-c 


we find 


(G a l ),r = ir\6(t - t')d r - t a (t,t'), 


(B.5) 

(B.6) 

(B.7) 

(B.8) 

(B.9) 

(B.10) 
(B. 11) 


where y cl = ri and y q = tq. 

The re-decoupled action ( |B .5| > is again quadratic in Majorana Grassmann variables tj a , enabling us to integrate 
them out and to obtain an effective action for the Q-fields. Thereafter, we can identify the saddle point and fluctuations 
of the effective Q-ac tion. 


iS[Q a ] = -Tr ( ^ log (G a {Q a )y l + -J dtdt’U clb (u / )M^Tv{rQ a (t,t'n b Qp(t'J)}. 


(B.12) 


Here Tr, denotes the trace in the Keldysh and time space and Tr is the trace in the Keldysh space. The saddle-point 
solution is found by expansion taking the linear order in 6Q. The solution must be stationary, depending only on the 
time difference: Gop(t, t') = Go pit - t’). We obtain the self-consistency equation 

Qo,a(t ~ t') = M a pGop(t - t’) . (B.13) 

In the high-temperature regime, T » TV, it is easy to obtain the self-consistent solution for the self-energy. In 
frequency space, one finds (summation over double indices assumed) 

Ka(“) = -\m„p f ^ (n K ( 0 J + Q)G A i8 (Q) + n s (w + O)G^(O)) 

= -2igT(l+0^ (B.14) 

= (oj- » (CJ + 2igT)- 1 , (B.15) 
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coinciding with the saddle-point result ( |35| > and the results of the perturbation theory ( |A.2[ i and ( ]A .3 [ >. 

To analyze the fluctuations 6Q, we expand the trace-log-term in ( |B.12| ) up to the second order in 6Q. With the 
help of the Fourier transform of 6Q, introduced as 

dQaico, v) = J dtdt'e^ (l+,,) e- i0J(t - n <5Q a (t,t'), (B.16) 

we rewrite the action in the form iS $q = iS ^ + iSfg. The first term originates in the expansion of the trace-log term 
of( |BT2l >: 


iS> 6Q 4/ dtdt ’ dt l dt l Tr {Go,q (t, t'l) s£ a (t', 0 } 

= \ I Tr {G 0 , a (fi - v/2) dt a (Q, v)G 0 , o (Q + v/2) 6± a (£l, -v)} . (B. 17) 

Here, for brevity, we expressed iS ^ via the self-energy fluctuations 

1 C* d 

v) = - - J ^Tifrtiin + cjyfsQafa-v)} . cb.is) 

The second term appeared after the Hubbard-Stratonovich transformation: 


m r 

= ~Y~ J dtdt’f\ ah (t, t')Tr {rSQ a (t, t’rfSQpit’, ?)} 

f ft ab (u>' - aj)Tr {y a 6Q a (o>, v)y h 6Q fj (oj', -v)} . (B. 19) 


K 
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In this Appendix we focus on the first term ( |B.17| >, which emerged from the expansion of the trace-log term. We 
give a detailed discussion of this term and conclude that it may be neglected as compared to the second term ( |B.19| >. 
This is done on the basis of a small-" expansion, justifying our course of action in Section 2.3 of the main text. In 
order to symmetrize and simplify the problem we parametrize fluctuations around the saddle point in terms of the 
mode R {1 \ symmetric with respect to time, and three antisymmetric modes R {2) ,R {3) and A" 41 : 


R ( a\t, f ) - R {2> (t',t), R (2) (t, t') - -R (2) (t\ t\ RV\t, n = -R%>(f, t), t') - -R^if, t) (B.21) 


?(Dr 


6Q a (u), v) 


? (2 V 


< 3) (tu,v) 
-iRa\u>, v ) ■ 
?(2)r 


■R ( a 4 \co,y) 

• R„\aJ, V ) 

?(3) 


iR ( J\oj, v ) ■ 


■ Ra 2) (o>, v) 


R^(aj,v)-R { f(oj,y) 


(B.20) 


( 3 ), 


?(4)r 




Clearly, (B. 19 1 is proportional to g essentially originating from the bath correlation function H. In (B. 17 1 , each c)X 


contains a factor of H, therefore the whole term appears to be of at least second order in g unless the Green’s functions 
yield an inverse factor g. The only combination of Green’s functions yielding 1 /g is G A (£1 + v/2)G K (£l - v/2) (or vice 
versa, R/A —» A/R). For example, one of the contributions to (|B. 1 1\ has the form 


f 


dcodu/'dvdQ 

&0* 


G a (G + v/2)G K (Q - v/2)n K (£2 + a>)n K (Q + oj')R d) (co, v)R' a> (oj\ -v), i, j = 3,4. (B.22) 




In this term, assuming » 4 gT to be constant (at low frequencies), the O integration of G R G A yields an inverse 
factor of g since 1/T = ( 2 gTY x . The structure of this term resembles the structure of the diagrammatic elements ( |A.7| ) 
discussed Appendix A However, in writing \B.22) we did not take into account the antisymmetry of R (3) and R i4> as 
explained in ( |B.21 1 . Due to the antisymmetry terms of the kind ( |B.22[ > cancel out. This is the same cancellation which 
we encountered in perturbation theory in Eq. ( |A.7| >, thus we conclude that the above mentioned divergent terms also 
cancel out in perturbation theory if symmetries are respected during the re-summation. 
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To substantiate our claim we provide here a rigorous analysis of ( |B.17| >. For this purpose we decompose c')X ( ,(Q, v), 
use the explicit form of W R ^ A and take advantage of the symmetry relations (|B.2 1 [) of R (l> . 


ds“(n,-v) = 


2n\ 2 

yi k (cj + Q) + n K (-a> + Q). 


if f (2f*n*JVv) 

/£( 

bX 2 '(Q, -v) = (Q,-v;<> -» -» -4 4) ) 

<® 22 (Q, -v) = (G, - v;^ -» -C;< 4) -> -4 4) ) 
In addition 


TP(w + Q) - n R (-to + G) 


(tf®(to,v)-^ 4 W)) 


i4 1) (tu,v) + 


n^(w + Q) - n K (-u + Q) 


R ( f\aj,v)-2gcjR^(cj,y) 


(% 


(B.23) 


n R/A (G + w) = ±g (G + to) - iD, n R (0 + to) = 2g(G + to) coth 


G + to 
2T 


(B.24) 


At low frequencies G «: T the terms containing the ’classical’ contribution I l A ~ 4 gT dominate over those containing 
n R/A . Indeed n R / A are important only in the ’quantum’ region S! > T, then II R / A / K oc gil. Here we recall the discussion 
about the RG in Section[2] We can disregard most of the ’quantum’ domain Q > T because these frequencies could 
be integrated out in the initial RG procedure. 

Where do large contributions to the O-integral in ( |B. 1 7[ ) come from? In the region Q ~ T Green’s functions 
generate a factor of 1/T 2 , the expression as a whole scales as g 2 and can therefore be neglected as compared with 


(B.19l, which scales as g. The remaining region to consider is Q ~ T <s T. There, terms containing the Keldysh 


Green function 


„ -2iV tanh 
G (n) = 


(B.25) 


get another order of g due to the hyperboli c tang ent in the small g expansion: tanh il/dT) ~ r/T = 2 g. Neglecting 
G A -terms we write the remaining terms of (B. 171 as 


G r bS 21 G R <ffi 21 + G r „ bS 22 G A v bS n + G a 6Z u GL bS 22 + G A ,M l2 G A WS 12 


(B.26) 


In the limit Q « T most prefactors in r)£ are linear in O. Considering the region Q ~ T = 2 gT the linear prefactor 
Q yields another order of g, and therefore the corresponding terms can also be neglected. Finally, the Q-indepcndcnt 
term aa>R {4) in <-)T n2 ' only appears combined with G R G R and G A G A , having both poles on the same side of the real 
axis. The integration by residue theorem yields zero. We conclude that all terms of ( |B. 1 7[ > are of higher order in g 
than those of ( B.19| >. This line of reasoning still holds if a magnetic field B is included in the problem provided that 
B <sz T. For larger fields, which are however still smaller than the temperature, the frequency in the hyperbolic tangent 
would essentially be replaced by B and thus yield a factor tanh B/(2T) ~ B/T > g. 

As a result we have confirmed the conclusion of the main text: it is justified to neglect the first term in the 
action ( |B.17| >, which was obtained from the expansion of the trace-log term around the saddle point. The action of 
dQ-fkictuations around the saddle point is governed by second term, ( |B.19| >, generated in the Hubbard-Stratonovich 
decoupling of the quartic Majorana term with Q-fields. 


Appendix C. Gauge freedom 

In this Appendix we explore the curious gauge freedom present in our problem. Interestingly, we find that the 
saddle-point solution ( |B.13| > and particularly the imaginary part of the self-energy acquire a gauge-field dependence. 
As a result, the physical Green function no longer coincides with the saddle point Green function. We find that 
fluctuations in turn become important, and their role is to compensate for the effect of the introduced gauge fields. 

We recapitulate the quartic term in the Majorana action (|B.3[), before the 0-fields were introduced. 


A M,int \Ho\ 


if 


dtdt'Tj IT 


ab 


( f > O A K a ri, a ) 


(C.l) 
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Due to the property rjj a = 0 of Grassmann variables, adding real/complex finite entries A x ,A y and A- on the diagonal 
of the matrix M does not change the action. The range of possible values of A a is limited by the constraints that M is 
invertible (or equivalently, det M + 0) and its eigenvalues have to be real. We interpret the A a as gauge fields, which 
may be fixed by some condition. 


M = 


ly -i 

1 A 


(C.2) 


The redecoupling of the quartic term with the help of complex bosonic fields Q„ via a Hubbard-Stratonovich transfor¬ 
mation is not modified, the action is still of the form given in the main text in Eq. ( |B.5| ). The only new feature consists 
in the diagonal non-zero entries of the matrix M. 


(C.3) 


iS [71a, Qa] = J dtdt' Urf a (t) (G a X T 1 h Jt’) - Tin ab (t, f) Q a a \t, f) M~pQ a p(t, f) 

Integrating out the Majorana fermions one again obtains ([32]) 

i W i p 

iS [Q a ] = -Tr f ^log (G a (Qa)T X + g J dtdt’U ah (t, l')M^Tr{y a Qa(t, , 


(C.4) 


and the saddle point equations 

Qo.aG, t’) = M a fjG()jj(t, t') 

(G 0 f)(t, t') - hiS(t - t')df - £(>,<*(*, t') 

Ka ^ f ') = -^ Tr {^n It, t')y b Qo,Jt', 0} • (C.5) 

We emphasize that this set of equations now depends on the gauge fields A a , introduced above, suggesting that there 
is not just one but instead a set of saddle points, characterized by the values of [A x ,A y ,A z }. This fact becomes obvious 
if we write down the self-consistent solution in the high-temperature regime explicitly: 


I R 


- -2 Ma P 


f ( n *( w + + n R (w + n)Gp(Q)) 



(C.6) 


As the imaginary part of now depends on the gauge fields, we can no longer identify it as a physically observable 
rate. For any finite [A x ,A y ,A z ] equation ( |42[ i fails and therefore the saddle point Green function and the physical 
Green function do not coincide, in contrast to the result ( |43| ) in the main text. In order to find the physical Green 
function, we have to reconsider fluctuations around the saddle points for finite gauge fields. 

To describe the fluctuations we evaluate the action (B.19), that is the leading second term of iSsQ■ To simplify the 
discussion we define 


< 1) (v)= f K ( J\r)= [ for i = 2,3,4. (C.7) 

Jo 2 n Jo In T 

Written in the matrix form in terms of (R^\, R ( ^\R \, 4) ), the leading terms in the high-temperature regime are 
O', j - 1 ^ 2, 3,4) 




(2 T 
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10/ \ 

r dv 

r° dojdcj' ,, 
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coco' 
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J 2n 

Jo art 
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3 T 
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(jJCl)' 

3 T 
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Joj 
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(0(1)' 

3 T ^ 


S r " ,(v) ( d X s a- v >. 


(C.8) 


(C.9) 
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0 

1 

3' 


(C.10) 


Thus, we obtain for the correlator of the fluctuations (r ( ^{v\ )R g> (v 2 )^ = D il j^ 2n6(vi + V 2 ). 

Knowing the propagator of fluctuations, we can compute the fluctuation-averaged Green function. As a starting 
point we use the self-consistent Dyson equation, which has to be satisfied by the correct Green function. 


(g~/ - X„) o G a = 1 . (C. 11) 

Before the average is actually performed causality does not necessarily apply, and we have to allow for a finite 
‘anti-Keldysh’ 22 component of the fluctuating self-energy. For the 12-component of the Green function, which will 
become retarded after the averaging, we obtain the following equation 


id,G 12 (t, f) = 6(t - t') ■ 


/ 


dt^ 2 a \t,h)G\ 2 {t u 


f) + 


I 


dhZ 22 {t,h)G 22 {h,t') 


(C.12) 


We have found that the propagator ( |C.10| ) does not depend on the frequency oj corresponding to the time difference 
t — t\. Therefore, we assume that X 21 (f, 1\) = 6(t - f |)X 21 (f) is local in time, but keep the dependence on total time for 
the fluctuation average later on. We also neglect X 22 G 22 , which is of higher order because X 22 ~ U K,A G R ^ A + 1 \ K G K is 
smaller than the ‘big’ combination Tl K G R ^ A . 


id,G 12 (t, t') = 6(t - f) + X 21 (f)G 12 (f, t’) 


(C.13) 


The above equation is solved by the following ansatz. The self-energy X 21 includes the constant saddle-point contri¬ 
bution and fluctuations. The saddle-point part is X^ and can be separated. The resulting prefactor is identified with 
the saddle-point Green’s function. 


- I2 (f,/') = -i®(t - f')exp j-z ^fiX 21 (fi)| = -i®(t - t')e n exp|-/J^ dt\dl, 21 (ti) 


= G R (t, t') exp 


-f 


dhSI^ih) 


(C.14) 


We treat the d()-fiuctuations using the usual Gaussian averaging procedure and use a ‘self-energy correlator’ for 
simplification. 


G>,0 = (G l2 (t,t')) 6Q = G^(f,f')|exp|-/ £ dt^l\h) 

t/f 2 t/f 3 (<5X 21 (f 2 )dX 21 (f3)) 


6Q 


= G o.J t ’ Oexp{-~ 


(C.15) 


The 21-component of the ‘self-energy correlator’ is decomposed into a bath correlator and fluctuating modes accord¬ 
ing to Eq. ( |B.18| ). 

<5X 21 (f,0 = - l -Tx{tKt,t')T X 8Q a {t' ,0} 

= K it - t) - \n K (t - t')R ( 3\f, t)- ] - (n*(f - f) + n A (t - o)4 3 V, t) 
- l -(n R (t-t')-n A (t-t'))R (A \t f ,t) 

=" K « - t')R<3\t',t) - ] - (n R u - f) - n \t - t'j)R%\t',t ) 
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X OO I 
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Ui 4 W 

Considering the exponent of (|C. 1 5|> we find 


(C.16) 


- £'dt 2 dt 3 (d^{ h )SJ} ct \ h )) 6Q = £'dt 2 dt 3 ^(sil , (y 2 )Sll , (v 3 )) 

(2 gT) 2 


-;v 2 f 2 -IV3?3 


-{t-t') (2i,0,0,-l) w £>; 


(0 

Q'Q' Q 


2i \ U) 


0 

0 

\-b 


= -gT{t-t')A a . (C.17) 


At the saddle point the imaginary part of the self-energy ( ]C.6| > was found to depend on the arbitrary constants 
A a . Hence it cannot correspond to the physical decay rate. We define the physical decay rates Iusing the physical, 
fluctuation-averaged 12-component of the Green’s function 

iG R a (t, t') = i (G^(t, t')) sQ = ©(t - t')e- r " ( - , - t ' ) (C. 18) 

and find 

r„ = - Im Z« a - gTA a = 2g T. (C. 19) 

In conclusion, we found that for arbitrary A a the physical decay rate is not given by the self-energy at the saddle 
point but rather by the decay rate V a of the fluctuation-averaged Green function, which is independent of A a . 

After having identified A a as kind of an arbitrary gauge, we can now choose A a = 0. Then, the saddle-point 
solution coincides with the correct decay rate, and corrections due to fluctuations cancel. In other words, the saddle 
point Majorana Green function coincides with the physical Green function. Analogously this applies to the self¬ 
energy. 
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